PROFESSIONAL  PAPER  361  /  September  1982 


METHODS  FOR 
GENERATING  AIRCRAFT 
TRAJECTORIES 


David  B.  Quanbeck 


dug 

/TH.  =  ■  ■  *'• 

(  "•'■■■ 1  &■ 

v  DEC  1  61982 


CENTER  FOR  NAVAL  ANALYSES 


This  documret  lias  boon  crp  moved  j 
for  pul'b-:  "1  Ua  I 

dUwrtoufK'n  i'  •.incited. 


82  12 


PROFESSIONAL  PAPER  361  /  September  1982 


METHODS  FOR 
GENERATING  AIRCRAFT 
TRAJECTORIES 


David  B.  Quanbeck 

riooem-?*-0' coa' 


Operations  Evaluation  Group 


CENTER  FOR  NAVAL  ANALYSES 


2000  North  Beauregard  Street,  Alexandria,  Virginia  22311 


ABSTRACT 


J. 

“^Methods  for  generating  three  dimensional  aircraft,  trajectories 
necessary  for  quantitatively  assessing  aircraft  tactics  are  documented 
in  this  report.  Elements  conventionally  used  in  modeling  aircraft 
motion  are  assembled  to  form  a  model  governing  aircraft  translation, 
fuel  use,  and  attitude.  Assumptions  on  the  functional  dependence  of 
the  aircraft  external  forces  and  specific  fuel  consumption  result  in  a 
system  of  seven  equations  ard  eleven  variables  governing  aircraft 
trajectories. 

To  provide  flexibility  in  prescribing  aircraft  trajectories,  the 
problem  of  solving  the  equations  is  formulated  for  five  separate  sets 
of  known  variables.  These  sets  include  variables  defining  aircraft 
controls,  velocity  attitude,  and  velocity  magnitude.  Extensions  to 
the  problem  formulations  allow  flight  path  normal  acceleration  to  be 
prescribed,  also.  A  method  to  prescribe  known  variables  is  presented 
that  ensures  continuous  aircraft  acceleration  and  angular  velocity. 
Numerical  integration,  finding  roots  of  equations,  and  interpolation 
of  function  values  are  required  to  solve  the  trajectory  generation 
problems.  Application  of  selected  algorithms  for  numerical  solution 
of  the  equations  is  discussed., 
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I.  INTRODUCTION 


The  mission  ef f Jctiveness  of  tactical  aircraft  can  be  assessed  by 
methods  ranging  from  using  mathematical  models  of  the  aircraft  and 
threat  weapon  systems  to  flight  testing  tactics  in  a  simulated  combat 
environment.  The  former  approach  is  useful  for  preliminary  assessment 
of  alternative  tactics  prior  to  employing  the  more  costly  flight 
testing  for  a  more  accurate  assessment.  Elements  used  to  mathemati¬ 
cally  model  engagements  between  an  aircraft  and  a  weapon  system  may 
include  a  model  that  governs  aircraft  motion,  models  of  aircraft 
subsystems  such  as  weapons  or  radar,  and  models  of  the  opposing  weapon 
systems  such  as  surface-to-air  missiles  and  radars.  Each  model 
reflects  the  inherent  capabilities  and  limitations  of  each  operational 
weapons  system  while  the  outcome  of  a  particular  engagement  is  an 
assessment  of  the  weapon  system's  overall  performance  given  the 
cactical  employment  of  the  systems  during  the  engagement. 

Aircraft  tactics,  in  particular,  vary  widely  due  to  the  varia¬ 
tions  in  the  trajectory  an  aircrew  can  rly  to  accomplish  a  mission  in 
addition  to  the  options  available  for  employing  any  of  the  subsystems. 
Quantitatively  describing  a  particular  aircraft  tactic  requires  the 
ability  to  generate  a  time  history  of  the  aircraft  trajectory  used 
during  the  tactic.  The  variables  describing  a  trajectory  that  are 
frequently  needed  in  assessing  aircraft  tactics  are  the  aircraft 
position,  velocity,  attitude,  and  fuel  use.  In  general,  a  model  for 
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generating  aircraft  trajectories  should  provide  for  aircraft  motion 
involving  arbitrary  three-dimensional  maneuvers  that  can  be  feasibly 
achieved  by  a  particular  aircraft  during  controlled  flight. 

This  report  documents  a  mathematical  model  and  solution  methods 
that  can  be  implemented  to  numerically  generate  aircraft  trajectories 
required  for  assessing  the  effectiveness  of  tactics.  The  model  is 
composed  of  elements  of  aircraft  dynamics  conventionally  used  in 
modeling  aircraft  trajectories.  Specifically,  six  scalar  equations  of 
motion  derived  from  the  vector  force  equation  expressing  Newton's 
second  law  and  an  equation  governing  aircraft  fuel  flow  comprise  the 
point-mass  model.  These  seven  equations  are  firBt  order  differential 
equations  governing  seven  variables  defining  the  aircraft's  velocity, 
position,  and  fuel  use.  Aerodynamic  forces,  engine  thrust  and 
specific  fuel  consumption  appear  in  the  equations.  This  information 
defines  the  inherent  capabilities  and  limitations  of  a  particular 
aircraft  in  terms  of  the  trajectories  it  can  feasibly  achieve.  Taking 
conventional  assumptions  for  the  functional  dependence  of  the  forces, 
four  additional  control  variables  defining  the  magnitude  and  orienta¬ 
tion  of  the  forces  are  needed  to  complete  the  set  of  variables  in  the 
model.  The  resulting  under-determined  system  of  seven  equations  and 
11  variables  can  be  solved  if  any  four  variables  are  prescribed  over 
time. 

Prescribing  four  of  the  variables  provides  control  over  the 
aircraft  motion  necessary  to  generate  a  particular  trajectory.  The 
choice  of  the  variables  that  are  prescribed  also  determines  the 
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procedures  required  to  solve  the  system  of  equations.  When  the  four 
control  variables  are  prescribed,  the  aircraft  motion  is  found  by 
numerically  integrating  the  equations  of  motion.  If  selected  state 
variables  are  prescribed,  unknown  control  variables  must  be  found  as 
roots  of  appropriate  governing  equations  with  the  remaining  equations 
integrated.  In  this  report,  the  solution  procedures  for  five 
different  sets  of  prescribed  variables  are  presented.  The  five  sets 
include  the  set  of  prescribed  control  variables,  two  sets  used  to 
prescribe  the  attitude  of  an  aircraft's  velocity  vector  with  selected 
controls,  and  two  sets  allowing  the  velocity  magnitude  and  attitude  to 
be  prescribed  with  a  selected  control  variable.  A  particular  set  of 
prescribed  variables  can  be  selected  according  to  which  set  allows  a 
particular  portion  of  a  trajectory  to  be  most  conveniently  defined.  A 
variation  of  the  problem  formulations  allows  the  flight  path  normal 
acceleration  to  be  prescribed  instead  of  one  of  the  angles  defining 
the  velocity  attitude. 

To  numerically  solve  the  equations,  the  variables  prescribed  as 
functions  of  time  can  be  constructed  using  arbitrary  functions  or  with 
a  procedure  presented  in  this  report  that  evaluates  the  variables  in 
terms  of  a  given  sequence  of  second  time  derivatives.  This  procedure 
ensures  that  linear  continuous  aircraft  acceleration  and  angular  velo¬ 
city  result  from  the  prescribed  variables.  Algorithms  for  integration 
and  finding  roots  of  equations  are  required  to  numerically  solve  the 
aquations  of  motion.  In  general,  each  derivative  evaluation  during 
numerical  integration  requires  finding  unknown  control  variables  as 


roots  of  algebraic  equations.  Additionally,  interpolation  between 
discrete  function  values  is  required  to  approximate  the  forces  and 
specific  fuel  consumption  appearing  in  the  equations. 

The  aircraft  trajectory  model  presented  in  this  report  includes 
elements  that  are  useful  for  investigating  other  problems  in  flight 
dynamics.  Flight  path  parameter  optimization  problem  formulations 
frequently  include  the  equations  of  motion.  Ibr  example,  a  parameter 
dependent  on  flight  path  variables  may  be  minimized  subject  to 
constraint  equations  which  include  the  equations  of  motion.  Adding 
the  moment  equations  governing  the  angular  motions  of  the  aircraft 
provides  a  model  that  may  be  used  to  investigate  aircraft  stability 
and  control  problems.  In  the  trajectory  model  presented  here,  the 
moment  equations  are  neglected  and  the  aircraft  angular  motions 
implied  during  a  trajectory  are  assumed  to  be  feasible.  The  moment 
equations  can  be  solved  given  the  angular  motions  during  a  particular 
trajectory  if  this  assumption  is  questioned. 

This  report  consists  of  four  sections  after  this  introduction. 

The  next  section  presents  the  model  governing  the  aircraft  motion. 
Included  in  this  section  are  definitions  of  variables  in  the  problem, 
development  of  the  scalar  equations  of  motion,  and  discussion  of  the 
assumptions  on  the  forces  and  specific  fuel  consumption.  The  third 
section  presents  the  individual  problems  formulated  with  the  different 
sets  of  prescribed  variables.  The  steps  required  to  solve  each 
problem  are  identified  for  both  the  general  problems  and  the  simpler 
cases  of  zero  sideslip  flight  and  symmetric  flight  in  the  vertical 
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plane.  Prescribing  the  velocity  attitude  in  terms  of  aircraft 
acceleration  is  considered  followed  by  a  discussion  of  limitations  and 
extensions  of  the  trajectory  model.  Numerical  methods  that  will  be 
used  to  implement  the  solution  of  the  equations  are  presented  in  the 
fourth  section.  These  include  a  method  for  prescribing  variables  and 
algorithms  chosen  for  integration,  root-finding,  and  interpolation. 
Special  consideration  is  given  to  application  of  the  algorithms  for 
solving  the  equations  in  the  trajectory  model.  The  last  section  of 
the  report  briefly  summarizes  the  model  and  presents  conclusions 
concerning  the  application  of  the  methods  for  generating  trajectories. 


II.  AIRCRAFT  MATHEMATICAL  MODEL 


This  section  presents  a  mathematical  model  governing  aircraft 
translation  and  fuel  use.  First,  three  reference  frames  are  intro¬ 
duced  and  the  transformations  between  the  reference  frame  coordinate 
axes  are  presented.  These  are  the  inertial,  the  wind,  and  the  body- 
fixed  reference  frames  used  for  representing  the  forces  acting  on  the 
aircraft  and  the  motion  of  the  aircraft.  Equations  for  calculating 
angular  velocities  of  the  moving  reference  frames  are  given  followed 
by  equations  for  calculating  the  aircraft  attitude.  The  scalar  force 
equations  of  motion  and  a  scalar  equation  governing  fuel  flow  are  then 
presented.  This  model,  selected  here  to  govern  aircraft  translation 
and  fuel  use,  is  the  point-mass  model  used  in  trajectory  analyses. 

The  model  neglects  the  equations  governing  the  aircraft's  angular 
motions  about  its  center  of  gravity.  This  section  ends  with  a 
discussion  of  the  trajectory  model.  Assumptions  are  taken  that  define 
the  dependence  of  the  force  and  specific  fuel  consumption  functions  on 
ctate  and  control  variables  in  the  model.  With  these  assumptions,  the 
trajectory  model  consists  of  seven  equations  and  11  variables. 

Defining  any  four  variables  as  known  functions  of  time  determines  a 
unique  trajectory.  Five  different  sets  of  known  variables  are  con¬ 
sidered  this  report  for  prescribing  aircraft  trajectories.  Discussion 
of  these  five  sets  and  the  corresponding  problem  formulations 
concludes  this  section. 


6 


This  section  includes  material  necessary  to  present  the  model 
governing  aircraft  trajectories  and  serves  to  document  the  equations 
comprising  the  model.  More  detailed  development  and  discussion  of  the 
elements  of  aircraft  dynamics  presented  here  can  be  found  in  [1]  and 
[21- 

reference  FRAMES 

Three  reference  frames  with  right-handed  coordinate  systems  will 
be  used  to  represent  aircraft  forces  and  motion.  These  reference 
frames,  the  inertial  frame,  the  wind  frame,  and  the  body  fixed  frame, 
are  described  below. 

inertial  frame,  Fj —  Newton's  laws  govern  the  motion  of  a  body 
with  respect  to  an  inertial  frame.  In  this  report,  the  inertial  coor¬ 
dinate  system,  xyz ,  is  assumed  fixed  on  a  flat  earth.  Acceleration 
of  the  aircraft  due  to  flight  over  the  rotating  curved  earth  is 
neglected  in  the  flat  earth  approximation.  Hie  increase  in  the 
acceleration  with  aircraft  speed  is  discussed  in  [1]  where  the  flat 
earth  approximation  is  considered  appropriate  for  flight  at  speeds 
below  about  Mach  3.  The  orientation  of  coordinate  axes  is  such  that 
z  is  positive  in  the  direction  of  positive  gravity,  g,  which  is 
also  assumed  constant.  The  directions  of  the  x  and  y  axes  are 
arbitrary . 

Wind  Frame,  Fw —  This  is  a  moving  frame  with  the  origin  of  the 
axes  xwywzw  fixed  at  the  aircraft  c.g.  The  axis  is  defined  to 

be  coincident  with  the  aircraft's  velocity,  V,  with  respect  to  the 


air  mass  (true  airspeed).  The  zw  axis  is  positive  in  the  lower  half 
of  the  aircraft  plane  of  symmetry,  downward  during  level  flight. 

Body-Fixed  Frame,  F^ — This  moving  frame  also  has  its  origin  at 
the  aircraft’s  c.g.  and  the  axes  xbybzb  are  with  respect  to 

the  aircraft.  The  axis  x^  is  defined  to  be  coincident  with  the 
zero-lift  longitudinal  axis  of  the  aircraft,  positive  forward.  The 
axis  z^  is  positive  in  the  lower  half  of  the  aircraft  plane  of 
symmetry. 

The  position  of  the  aircraft  c.g.  in  the  inertial  frame  will  be 

noted  by  the  vector  -  [x^,  y-j-,  Zjjfc.  The  aircraft’s  inertial 

"  •  •  •  *. 

velocity,  the  time  derivative  of  Xj,  is  then  Vj  *  [xi ,yj * zi ]  » 

If  air  mass  has  a  velocity  W  with  respect  to  the  inertial  frame, 
then  Vj  *  V  +  W.  In  the  following  development,  W  will  be  assumed 
constant  in  time  and  space. 

The  coordinate  axes  of  the  moving  reference  frames  are  displaced 
by  translation  and  rotation  from  the  inertial  axes.  Components  of  the 
same  vector  observed  from  two  parallel  coordinate  systems  are  equal, 
but  angular  orientations  of  the  moving  coordinate  systems  need  to  be 
defined  to  develop  matrices  for  transforming  vectors  between  rotated 
coordinate  systems.  The  orientation  of  the  vind  and  body-fixed  axes 
are  shown  in  figure  1  where  x’  y*  z*  is  a  set  of  axes  parallel  tc 
the  inertial  frame. 

Three  Euler  angles  designate  the  orientation  of  thr  wind  cues 
from  the  inertial  axes.  HL>  ;t:s  x*  and  y'  are  rotated  about  the  z' 
axis  to  form  an  intermediate  coordinate  system  x^,y^,z’  where  x^ 


is  coincident  with  the  projection  of  V  in  the  horizontal  plane.  Th 
angle  of  rotation  is  the  velocity  yaw  angle,  \p.  A  rotation  0  ,  the 
velocity  pitch  angle,  about  y^  carries  x^  to  xw,  coincident 
with  V,  resulting  in  a  second  set  of  intermediate  axes  xw>yi>zi* 
The  velocity  roll  angle,  <f>  ,  is  the  final  rotation  about  xw  to 
carry  z 1  into  the  aircraft  plane  of  symmetry  thus  forming  the  wind 
axes  x  y  z  .  Each  of  the  Euler  rotations  is  a  rotation  of  two  axes 
in  a  plane,  so  three  vector  coordinate  transformations  about  a  single 
axis  occur  in  sequence.  Theso  rotations  result  in  a  matrix,  L^,  to 
transform  a  vector  expressed  on  the  inertial  axes  to  the  same 

vector  Aw  **  L^Aj  expi  eased  in  Fw* 


COS0COSlp 


cosOsirap 


-sin0 


L 


wl 


sin<psin0cosip  sin<psin9sintp 
-cosips  imp  +cos<pcos\p 


sirwpcos0 


(1) 


cos<psin0cosip 
+s  i  n<p  i  i  nip 


c^4  sin9  simp 
•  sinpr.osip 


COSlj)COS0 


LwI  is  an  orthogonal  matrix,  its  inverse  is  equal  to  its  transpose. 
Thus,  the  transformation  from  Fw  to  Fj  is  gi'-en  by  A^  *  ^iw^w» 
where 

The  orientation  of  the  body-fixed  axes  with  respect  to  the  wind 
axes  is  also  defined  by  two  Euler  angles.  A  rotation,  -0  ,  about 


z  results  in  the  x2^hzw  intermediate  axes.  The  quantity  0  is 
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the  sideslip  angle  and  the  *2ZW  P^ane  *s  also  :he  aircraft  plane  of 
symmetry.  Rotating  the  coordinates  in  this  plane  about  y^  through 
the  angle  of  attack  a  yields  the  aircraft  body-fixed  axes 
The  orthogonal  transformation  matrix  resulting  from  tnese  two 
rotations  is: 


cosctcosg 

sing 

sinacosg 


-coscxsing 

cosg 

-sinasing 


-sina 

0 

cosa 


(2) 


Expressions  for  the  angular  velocities  of  the  moving  coordinate 
systems  are  developed  next.  Of  particular  interest  is  the  angular 
velocity  of  the  axes  x  y^z^  to  be  used  to  express  the  aircraft 
inertial  acceleration  in  this  coordinate  system.  Hie  aircraft  angular 
velocity,  with  components  being  the  aircraft  yaw,  pitch,  and  roll 
rates,  is  the  angular  velocity  of  xbybzb*  Monitoring  the  value  of 
.hese  components  will  be  useful,  especially  during  maneuvers  involving 
rapid  changes  of  aircraft  attitude. 

The  total  angular  velocity  vector  of  a  moving  coordinate  system 
with  respect  to  the  inertial  frame  is  the  sum  of  the  angular  velocity 
vectors  due  to  the  time  rate  of  change  of  each  of  the  Euler  rotations. 

For  the  axes  the  angular  velocity,  u>wt ,  due  to  the  rates 

•  * 

0  and  is 


(3) 
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Here,  and  k'  are  the  unit  vectors  on  the  y^  and  z'  axes, 

respectively.  Observing  that  k'  -  -sin9i  +  cosQk, ,  where 

w  l  w 

and  are  unit  vectors  on  and  y^,  gives  w  »  expressed  on 

Wi  as 


w' 


s  i  nQ 
9 

^cos0 


(4) 


Calculating  the  angular  velocity  of  the  body  axes  requires  adding 

•  •  • 

the  components  due  to  4> ,  —6 ,  and  a  to  the  two  components  8  .mined  in 
equation  3.  FLrst,  the  angular  velocity  of  the  wind  axes  xwywzw  can 
be  found  with  the  components  expressed  in  terms  of  the  wind  axes 
coordinates.  Employing  the  approach  used  to  obtain  equation  4,  the 
resulting  wind  axes  angular  velocity  is 


U) 

w 


4>  — ^s  I  nQ 

•  • 

0co8<j>  +  <J>sin$cosQ 

•  • 

-0sin4>  +  i}»cas$cos8 


(5) 


•  • 

The  angular  velocity  vectors  due  to  -0  and  a  are  next 
expressed  in  terms  of  the  body-fixed  coordinate  system.  Summing  this 
vector  with  u)^  transformed  to  the  body-fixed  frame  gives  the 
following  expression  for  calculating  the  aircraft  angular  velocity. 


V 

' 

• 

3slnot 

“b  3 

qb 

KX 

• 

a 

• 

rb 

-Scosa 

(6) 
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The  components  pb,  qfe,  and  rfe  are  the  aircraft  yaw,  pitch  and  roll 
rates,  respectively. 

Knowledge  of  an  aircraft's  attitude  with  respect  to  the  inertial 
axe«  is  frequently  important  in  assessing  a  trajectory.  An  aircrew's 
f ield-of-view,  aircraft  sensor  and  weapons  employment  envelopes,  and 
the  aircraft's  aspect  as  seen  by  a  second  observer  are  examples  of 
quantities  dependent  on  the  aircraft  attitude.  The  orientation  of  the 
body-fixed  axes  from  the  inertial  axes  can  be  defined  by  three  Euler 
angles  0^,  and  <J>b  which  are  the  aircraft  yaw,  pitch,  and  roll 

angles.  These  angles  are  analogous  to  the  wind  axes  Euler  angles  and, 
therefore,  a  transformation  matrix,  L^j,  results  identical  to 
equation  1  except  that  @b,  and  4>b  replace  the  corresponding  wind 

axes  Euler  angles.  To  calculate  the  aircraft  attitude,  terms  of 
can  be  equated  to  terras  in  the  matrix,  {*^}  »  equal  to  the  product 
Lbw*LwI‘  This  results  in  the  expressions  below  for  the  body-fixed 
Euler  angles. 


-1  *12 
4V  "  tan  - 

b  *11 


(7a) 


-sin'1*,,,  -it/2  <  6U<  n/2 

1 J  D 


(7b) 


,  ,  -1  23 

♦  K  “  tan  j — 
b  *33 


(7c) 


Above,  the  terras  of  l^^jl  are  the  results  of  evaluating  trigono¬ 
metric  functions  of  a,  8»  and  the  wind  axes  Euler  angles.  The  signs 


of  the  arguments  In  the  Inverse  tangent  functions  above  will  determine 
the  appropriate  quadrants  of  'J’b  and  ♦t* 

EQUATIONS  OF  MOTION 

This  section  presents  the  scalar  equations  of  motion  governing 
aircraft  translation.  First,  three  scalar  force  equations  governing 
aircraft  velocity  are  formulated  on  the  moving  axes, 
aircraft  position  In  the  inertial  frame  is  governed  by  three 
additional  equations.  A  single  equation  governing  aircraft  fuel  flow 
completes  the  model. 

A  derivation  of  the  force  equation  governing  aircraft  motion  is 
presented  in  [2].  The  resulting  vector  equation  consistent  with  the 
flat  earth  approximation  is 

T  +  A  +  mg  «  maj.  .  (8) 

where  T  -  aircraft  thrust 

A  *  aerodynamic  force 
g  ■  acceleration  of  gravity 
m  ■  aircraft  mass 

■  inertial  acceleration  of  the  ai. 'craft  mass  center. 

•  •  • 

ihree  scalar  equations  for  the  time  derivatives  V,  i|»,  and  0  will 
be  found  next  by  expressing  the  components  of  the  equation  8  on  the 
moving  coordinate  system  xwyizi*  *n  doing  so,  each  of  three 


derivative  terms  will  appear  in  only  one  equation,  a  convenient  form 


for  numerical  integration.  However,  the  components  of  T  and  A 
will  be  summed  first  on  the  wind  axes  and  then  transformed  to  x^j^. 

Net  thrust,  T,  is  assumed  to  act  in  the  aircraft  plane  of 
symmetry,  x^z^,  at  a  angle  e  elevated  from  the  aircraft 

longitudinal  axis,  x^.  Transforming  the  thrust  v  ctor  from  the  body 
fixed  axes  to  the  wind  axes  gives 


Tcose 

-  *i 

TcosBcos(ct+e ) 

T  -  L  , 
wb 

0 

m 

-Ts  in0  cos(a+e  ) 

-Taint 

-Tsin(a+e ) 

The  three  components  of  T  above  will  be  noted  as  T  ,  T  ,  and 

T  .  The  aerodynamic  force  vector  components  are  defined  in  the  wind 
w 

axes  as  drag,  side  force  and  lift  on  the  xw,  yw„  and  zw  axes, 
respectively.  All  three  components  are  assumed  to  act  in  negative 
direction  of  their  respective  axes  giving 


(10) 


The  wind  axes  are  obtained  from  by  a  single  rotation  4> 

the  xw  axis.  Then  T+A  is  transformed  to  x^^z^  by 


"l 

• 

o 

c 

T  -D 

X 

w 

T  -D 

X 

w 

T  +  A  = 

0 

cos  —si  n4> 

T  -C 
*w 

m 

cos<fr(T  -C)  -  sin<|>  (T  -L) 
y  z 

w  w 

0 

s  i  n«J>  C0  5(|i 

T  -L 
z 

w 

simJ»(T  -C)  +  cosi>(T  -  L) 
^w  Zw 

about 


(ID 


The  aircraft's  weight,  mg,  has  only  one  non-zero  component,  expressed 
on  the  coordinates  x'y'z',  equal  to  rag£' .  As  in  the  development  of 
equation  4,  lc'  can  be  replaced  by  -sin9 i^+cosOk^  .  Therefore, 
aircraft's  weight  expressed  on  xw>izi  is  given  by 


mg 


mg 


-sin0 

0 

COS0 


(12) 


Finding  the  aircraft's  inertial  acceleration  remains  to  complete 
the  scalar  force  equations.  Recalling  ■  V  +  W  and  that  W  is 
assumed  to  be  constant,  the  acceleration  is  then  the  time  derivative 
of  V,  dV/dt,  observed  from  the  inertial  frame.  However,  the 
acceleration  vector  is  to  be  expressed  on  the  moving  frame  wy^z^. 

In  terms  of  V  and  ooTV «  as  observed  from  xwyizl»  the  inertial 
acceleration  expressed  on  x^^z^  is  the  sura  of  the  derivative  of  V 
and  ui^xV.  By  definition,  V  has  one  component  ou  Vlw» 

Using  equation  4  defining  the  components  of  ww»  gives  the  inertial 
acceleration  vector  on  xwy\zi  as 
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(13) 


Having  developed  the  components  of  the  vectors  in  equation  5 


expressed  on  xwyizi*  three  equivalent  scalar  equations  can  be 
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written.  Solving  for  the  derivative  terms  in  the  acceleration 
components  gives 


V  -  -  (T  -  D]  -gsinQ 

m  x 

w 


(Ua)  ; 


i 


*  -  ;sibir  'cos*(Ty  - c) '  s1i4(t2  _  L)1 

w  w 

0  -  [sin4>(T  -  C)  4-  coaiji  (T^  -  L)  ]  -  &  cos0  (14c) 

m  y  w  w 


These  three  first  order  ordinary  differential  equations  govern  the 


magnitude  and  orientation  of  V  at  each  time  instant.  The  equations 
governing  the  aircraft's  position  in  the  inertial  frame, 


XI  "  txI» 

yI’  ziJt» 

are  obtained  from  dX^/dt  * 
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coordinate 

system, 

Vx  is  the  sum  of  and 

the 

velocity  of  the 

air  mass. 

W  -  [wx. 

Wy,  Wz]t. 

• 

XI 

■  VcosGcos^  +  w^ 

(15a) 

• 

yI 
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ZI 

=■  -VsinO 

+  w 
z 

(15c) 

To  relate  fuel  use  to  aircraft  motion,  one  additional  equation 
will  complete  the  model  of  aircraft  ~lght  expressed  by  equations  14a 
through  15c.  As  fuel  is  burned  to  produce  thrust,  the  mass  of  the 
aircraft  decreases.  A  discussion  of  turbojet  and  turbofan  engines  is 


presented  in  [2]  giving  the  following  general  representation  for  the 
mass  flow  thrust  relationship. 

m  =■  -cT/g 

(16) 

c  =  specific  fuel  consumption 

DISCUSSION  OF  THE  EQUATIONS  OF  MOTION 

Equations  14a  through  16  constitute  the  model  of  aircraft  motion 
assumed  for  the  rs-iw-iinder  of  this  report.  Seven  first  order  ordinary 
differencial  equations  govern  the  seven  state  variables  (V,  i|t ,  0,  Xj, 
yI,  Zj,  in).  The  remaining  variables  iuc’ude  the  angles  (a,  0,  <{>), 
the  thrust  and  aerodynamic  forces,  and  the  specific  fuel  consumption. 
In  the  discussion  below,  assumptions  are  made  about  the  functions 
defining  the  forces  and  the  specific  fuel  consumption.  These 
assumptions  reduce  the  number  of  variables  in  the  model.  Then,  five 
alternative  problem  formulations  are  presented.  The  alternatives 
arise  from  prescribing  different  sets  of  known  variables  over  time  and 
solving  for  the  remaining  variables. 

Aerodynamic  forces  are  usually  represented  by  dimensionless 
coefficients  obtained  by  dividing  the  forces  by  the  product  of  dynamic 
pressure  and  a  reference  area.  The  coefficients  are  functions  of 
aircraft  attitude,  shape,  Mach  number,  and  Reynold's  number.  Varia¬ 
tions  in  the  forces  due  to  Reynold's  number  effects  usually  can  be 
neglected,  see  [1].  Furthermore,  the  drag  and  lift  will  be  assumed 
independent  of  the  sideslip  angle  and  the  side  force  will  be  assumed 
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independent  of  the  angle  of  attack.  With  these  assumptions  the  drag, 
side  force,  and  lift  coefficient  functions  are 


Cd(K,  a) 


2D/(p72S) 


CC(M,  (3)  -  2C/(pV*S) 


Cl(M,  a)  -  2L/(pV  S) 


(17a) 


(17b) 


(17c) 


where 


M 


a 


P 

S 


V/a  =»  Mach  number 
speed  of  sound 
atmospheric  density 
reference  area  (wing  area) 


The  three  force  coefficient  functions  will  be  numerically  approximated 
by  interpolating  between  values  of  the  functions  stored  in  two  dimen¬ 
sional  data  sets.  This  data  will  apply  to  a  given  aircraft  shape. 
Variation  in  aircraft  shape,  by  carrying  external  stores  or  changing 
wing  sweep  angle,  may  significantly  change  a  force  coefficient 
function.  Either  corrections  to  the  function  values  or  increasing  the 
dimension  of  the  domain  of  the  function  are  required.  Extensions  of 
the  functions  can  be  determined  f  r  .•  particular  aircraft  and  will  not 
be  considered  further  here. 

Thrust  and  specific  fuel  consumption  of  turbo' an  and  turbojet 
engines  can  also  be  represented  by  dimensionless  coefficients 


presented  in  [2].  The  coefficients  are  assumed  independent  of 
Reynold's  number  and  the  angles  a  and  0.  Engine  performance  is 
then  defined  by  the  thrust  and  specific  fuel  consumption  coefficients 
below. 

VM-  nc)  ■  w  <18‘> 

e 

ca* 

Kc(M’  V  -  IT  <18b> 

where  p  *  atmospheric  pressure 
S„  *  reference  area 

a*  -  speed  of  sound  at  the  tropopause 
nc  ■  corrected  engine  speed  ■  na*/(anmax) 
n  *  engine  rotor  rpm 

Like,  the  aerodynamic  force  coefficients,  the  above  functions  can  be 
approximated  by  two  dimensional  data  sets  for  a  given  aircraft.  Also, 
the  functions  can  be  used  in  an  arbitrary  atmosphere  since  the  speed 
of  sound  is  proportional  to  the  square  root  of  the  temperature  of  the 
air.  However,  in  [3],  a  discussion  of  engine  modeling  states  that 
approximating  T  and  c  to  be  strictly  proportional  to  p  and  a 
will  result  in  errors.  These  errors  result  from  Reynold's  number 
effects  and  deviations  from  the  assumed  thermodynamic  properties  of 
the  engine  airflow.  If  the  desired  accuracy  requires  these  effects  to 


be  included,  then  pressure  altitude,  with  a  standard  day  atmosphere 


assumed,  can  be.  added  to  the  domain  of  the  functions.  If  so, 
equations  18a  and  18b  can  still  be  used  to  approximate  engine 
performance  in  atmospheric  conditions  deviating  from  a  standard  day  at 
l  given  pressure  altitude.  In  this  report,  equations  18a  and  18b  will 
be  assumed  to  model  aircraft  engine  thrust  and  specific  fuel 
consumption. 

The  aircraft  model  defined  by  equations  14a  through  16  can  now  be 
interpreted  in  view  of  the  assumed  functional  dependence  of  the  forces 
and  specific  fuel  consumption.  Assuming  the  properties  of  the 
atmosphere  are  known  functions  of  altitude,  -zx>  the  m°del  consists 
of  seven  equations  determined  by  eleven  variables.  The  variables  can 
be  grouped  into  the  seven  state  variables  governed  by  the  seven 
differential  equations,  (V,  ip,  6,  Xj,  yIf  Zj,  m) ,  and  four  control 
variables  (a,  8,  <|> ,  nc).  In  an  aircraft  dynamics  model  including  the 
moment  equation,  as  in  [1],  a,  0,  and  <j>  appear  as  state  variables 
primarily  controlled  by  deflections  of  the  elevator,  rudder  and 
ailerons,  respectively .  The  corrected  engii  a  speed  is  controlled  by 
the  throttle  position. 

With  seven  equations  and  eleven  variables,  the  problem  of 
generating  an  aircraft  trajectory  can  be  solved  numerically  by  pres¬ 
cribing  any  four  variables  as  known  functions  over  time  and  defining 
initial  conditions  for  the  state  variables.  Of  all  the  possible 
combinations  of  four  known  variables,  five  different  combinations  will 
be  considered  in  this  report.  Hie  first  set  will  be  defined  by 
prescribing  the  four  control  variables  (a,  0,  ,  nQ).  Then  the  seven 


equations  are  integrate'  _o  find  aircraft  velocity,  position  and  fuel 
use.  Two  more  combinations  of  known  variables  considered  are  defined 
by  prescribing  the  attitude  of  the  velocity  vector,  the  engine 
control,  and  either  <j)  or  3*  Numerical,  solutions  given  either  of 
these  sets,  (^ ,  9,  <J> ,  a  )  or  (tj» ,  9,  3,  nc;,  requires  solving  equations 
14b  and  14c  for  the  roots  (ot ,  3''  or  (a,  <f>  >  and  integrating  the 
remaining  equations.  The  last  two  sets  of  known  variaKles  are  defined 
by  prescribing  the  velocity  vector  and  either  <j>  or  3,  that  is 
(V,  K  0,  <>)  or  (V,  ,  0,  3)*  Then  the  three  force  equations,  i4a, 

14b,  and  14c,  are  solved  for  the  roots  (a ,  3 ,  nc)  cr  (a,  <|>,  nc) . 

Inte'  ating  equations  15a  through  16  gives  aircraft  positlo  1  and  fuel 
use. 

In  summary,  any  one  of  v  e  sets  of  variables  can  be  prescribed 
to  formulate  the  problem  of  solving  aquations  14a  through  16.  .'I though 
any  one  set  is  sufficient  to  solve  for  any  feasible  aircraft 
trajectory,  these  alternative  problem  formulation-  allow  flexibility 
in  prescribing  aircraft  motion.  With  this  approach,  the  sever, 
equations  can  be  solved  for  a  series  of  aircraft  maneuver t  comprising 
a  complete  trajectory.  Each  maneuver  can  be  generated  by  prescribing 
the  set  of  variables  vith  which  the  maneuver  is  most  conveniently 
defined.  Example  flight  cot  ditions  easily  defined  with  the  different 
sets  of  known  variables  are  given  in  the  next  section. 


III.  TRAJECTORY  GENERATION  PROBLEM  FORMULATIONS 


Etch  set  of  described  variables  discussed  at  the  end  of  the  last 
section  determines  a  different  formulation  for  the  problem  of  numeri¬ 
cally  solving  equations  14a  through  16.  The  five  different  problems 
are  discussed  in  this  section  in  terms  of  the  procedures  necessary  to 
solve  the  problems.  When  state  variables  are  known  variables,  control 
variables  must  be  found  in  general  by  solving  simultaneous  equations. 
Fbr  zero  sideslip  flight  and  symmetric  flight  in  the  vertical  plane, 
particular  cases  of  interest,  solving  the  sets  of  algebraic  equations 
simplifies.  In  addition  to  discussing  the  different  problem 
formulations,  examples  of  specific  aircraft  maneuvers  conveniently 
formulated  in  each  case  are  given.  A  trajectory  can  be  constructed  by 
a  sequence  of  maneuvers  each  prescribed  with  one  set  of  known 
variables.  Later  in  this  section,  expressions  for  prescribing 
velocity  attitude  angles  in  terms  of  the  acceleration  normal  to  the 
flight  path  are  presented.  Final  'y,  limitations  of  the  aircraft 
trajectory  model  are  considered  followed  by  a  brief  discussion  of 
aircraft  dynamics  problems  that  may  be  solved  by  alternative 
formulations  or  extensions  of  the  model. 

PRESCRIBED  CONTROLS 

In  this  case,  the  control  variables  (a,  3,  <j>  ,  nc)  are  defined 
as  known  functions  of  time.  Given  initial  values  for  the  state 
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variables,  equations  14a  through  16  are  integrated  over  time  giving 
aircraft  velocity,  position  and  fuel  use.  To  evaluate  the  derivatives 
of  the  state  variables  for  numerical  integration  the  variables  and 
functions  appearing  in  equations  14a-c  and  16  must  be  evaluated. 

Several  specific  steps  are  required  prior  to  calculating  the 
state  variable  derivatives  at  any  time,  t.  First  the  known 
variables,  ot(t),  B(t),  4>(0»  and  n£(t)  have  to  be  evaluated.  For 
example,  they  may  be  numerically  evaluated  from  analytic  functions 
provided  for  a  particular  problem  or  evaluated  by  interpolation 
between  G^screte  points.  Then  the  properties  of  air  at  the  current 
altitude  and  the  Mach  number  must  be  calculated.  The  aerodynamic 
force,  thrust,  and  specific  fuel  consumption  coefficients  are 
evaluated  by  interpolation  between  stored  discrete  function  values. 
These  steps  will  be  required  in  all  five  problem  formulations. 

Finally,  equations  14a  through  16  for  the  state  variable  derivatives 
are  evaluated.  The  numerical  algorithms  to  accomplish  each  step  are 
presented  later  in  a  separate  section. 

When  3*0  the  side  force  is  also  zero  and  this  condition 
together  with  nin<J>  ■  0  results  in  symmetric  flight  in  the  vertical 
plane.  Equation'.  14b  (and  15b  if  ^  ■  0,  Wy  *  0)  are  removed  from 
the  system  of  equations  reducing  the  integration  problem.  These 
assumptions  will  be  convenient  in  calculating  fuel  use  over  two- 
dimensional  trajectories  with  (a,  nc)  prescribed,  for  example, 
simply  as  constants.  For  flight  in  three  dimensions,  control 
schedules  may  be  formulated,  for  example,  by  modifying  control 
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schedules  found  from  solving  the  equations  with  a  different  set  of 
known  variables.  A  general  extension  of  the  problem  formulation  would 
be  to  evaluate  controls  during  the  integration  based  on  the  state  of 
the  aircraft  motion.  Ebr  example,  fuel  use  during  constant  lift 
coefficient  trajectories  can  be  evaluated  by  appropriate  selection 
of  a  as  the  integration  proceeds. 

PRESCRIBED  VELOCITY  ATTITUDE  AND  ENGINE  CONTROL 

Selecting  either  of  the  sets  (i|i,  B,  n£l  3)  or  (<J>,  8,  n£,  <f>)  as 

the  known  variables  will,  in  general,  require  solving  equations  14b 

and  14c  for  the  unknown  controls  a  and  f  or  p.  Initial  values  of 

the  remaining  variables  (V,  xj,  Zj,  yj,  m)  must  also  be  defined.  At 

any  time  t  the  first  step  in  the  solution  is  to  evaluate  i|»(t),  0(t), 

nc(t),  and  B(t)  or  $(t).  Values  for  the  derivatives  of  the  known 
•  • 

state  variables  <Kt)  and  9(t)  must  also  be  evaluated.  Then,  the 
unknown  controls,  (a,  $)  or  (a ,  0 )  that  satisfy  equations  14b  and 
14c  have  to  be  found.  The  Newton-Raphson  root-finding  algorithm 
selected  for  this  step  in  the  solution  is  presented  later.  In  general, 
the  equations  will  have  to  be  solved  simultaneously,  and  repeated 
evaluations  of  the  forces  in  the  equations  are  required  during  the 
iterative  root-finding  algorithm.  Once  the  roots  are  found,  the 
derivative  of  the  state  variables  in  equations  14a  and  15a  through  16 
are  evaluated  for  numerical  integration. 

Two  special  cases  of  interest  are  symmetric  flight  in  the 
vertical  plane  and  zero  sideslip  flight.  Symmetric  flight  in  the 


vertical  plane  results  when  the  set  of  known  variables  includes 

• 

■  0  and  sin$  -  0  or  0  ■  0.  Then,  the  problem  reduces  to 
finding  a  as  the  root  of  14c  and  Integrating  14a  and  15a  through 
16.  Defining  and  w^  equal  to  zero  further  simplifies  the 
problem  deleting  15b  from  the  integrated  equations.  Zero  sideslip 
flight  results  from  defining  0*0  and  the  solution  simplifies  as 
<j>  can  be  solved  independently  of  a.  Solving  equation  14b  and  14c 
for  4>  gives 

4>  ■  tan  *[tcos0/(0  +^-cos0)]  (19) 

The  appropriate  quadrant  of  $  is  determined  by  the  signs  of  the 
numerator  and  denominator  of  the  arctangent  argument.  Equation  14c  is 
then  solved  for  the  root  a  and  integration  of  14a  and  15a  through  16 
proceeds  as  in  the  general  case. 

Selecting  either  (<Ji,  9,  nc,  0)  or  ,  0  ,  nc,  <j>)  as  known 
variables  will  be  useful  for  finding  aircraft  velocity,  position,  and 
fuel  use  during  flight  easily  described  by  the  direction  of  the  air¬ 
craft  velocity  vector.  Symmetric  flight  in  the  vertical  plane  can  be 
assumed  when  generating  aircraft  maneuvers  such  as  level  acceleration, 
climbs,  and  pull-ups  from  level  flight.  Zero  sideslip  flight  may  be 
assumed  in  modeling  turning  flight  as  in  level,  climbing,  or 
descending  turns.  Prescribing  non-zero  0  can  generate  motion  such 
as  turns  with  sideslip.  Prescribing  $  would  be  useful,  for  example, 


when  rolling  the  aircraft  to  the  Inverted  attitude  prior  to  the 
transition  from  a  climb  into  a  dive. 

PRESCRIBED  VELOCITY  VECTOR 

Prescribing  the  magnitude  and  direction  of  the  aircraft's 
velocity  vector  over  time  is  possible  when  (V,  9,  \p ,  3)  or  (V,  0, 
iJj  ,  ij> )  are  selected  as  known  variables.  Given  initial  values  for  the 
remaining  state  variables,  the  solution  requires  solving  14a,  14b,  and 
14c  for  the  roots  a,  n  and  4>  or  0  and  then  integrating  equations 
15a  through  16.  The  solution  steps  in  this  case  differ  from  those  of 
Luo  previous  problem  formulation  in  that  three  simultaneous  equations, 
instead  of  two  equations,  must  be  solved  followed  by  integration  of 
15a  through  16.  The  conditions  determining  symmetric  flight  in  the 
vertical  plane,  and  zero  sideslip  flight  in  the  preceding  case  also 
apply  in  these  problem  formulations.  Ebr  both  symmetric  and  zero 
sideslip  flight,  the  root-finding  problem  reduces  to  solving  14a  and 
14c  simultaneously  for  a  and  nc .  Fbr  zero  sideslip  flight, 
equation  19  gives  the  value  of  4>  needed  to  solve  14a  and  14c. 

Prescribing  the  velocity  vector  in  particularly  useful  in  solving 
for  the  controls  and  fuel  use  associated  with  steady  flight  conditions 
(V  *  0)  which  are  conveniently  expressed  by  (V,  9  ,  ij» ,  0 )  or 
(V,  0,  ) .  Constant  velocity  turns,  climbs  and  level  flight  are 

typical  examples.  Specifying  the  velocity  vector  also  provides  a 
convenient  way  to  start  a  trajectory  from  a  steady  flight  condition 
prior  to  generating  aircraft  maneuvers. 
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The  .  ive  different  sets  of  known  variables  considered  have  been 
discussed  in  terms  of  three  general  problem  formulations.  These  three 
problems  are  characterized  by  particular  equations  that  must  be  solved 
for  unknown  control  variables  and  those  to  be  Integrated  to  find  the 
state  variables.  Methods  for  prescribing  the  known  variables  as 
functions  of  time  have  not  been  discussed.  A  method  useful  for 
prescribing  any  of  the  known  variables  is  presented  in  the  next 
section.  An  extension  to  the  general  problem  formulations  resulting 
when  0  and  are  prescribed  variables  will  be  presented  next  that 
will  allow  additional  flexibility  In  prescribing  aircraft  maneuvers. 
The  extended  formulations  allow  6  or  |  to  be  replaced  by  the 
flight  path  normal  acceleration  in  the  sets  of  prescribed  variables. 

FLIGHT  PATH  NORMAL  ACCELERATION 

Aircraft  motion  involving  changes  in  velocity  yaw  or  pitch  angles 
requires  forcea  acting  on  the  aircraft  often  significantly  larger  than 
the  forces  encountered  during  steady  flight.  These  forces  should  not 
exceed  aircraft  structural  limits  and  the  resulting  acceleration 
should  not  exceed  acceleration  tolerable  by  aircrews.  The  magnitude 
of  the  acceleration  normal  to  the  aircraft  flight  path  (in  the  y^z^ 
plane)  can  be  found  from  the  acceleration  components  on  y^  and 
as 


„  *2  ,  • 2  . 2  xl/2 

■  V(9  +  co 80  ) 


a 


(20) 
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Prescribing  a(t)  in  place  of  either  t|>(t)  or  9(t)  is  a  convenient 
way  to  characterize  certain  maneuvers,  especially  those  to  be  limited 
by  acceleration  magnitude.  Depending  on  whether  0(t)  or  <|>(t)  is 
prescribed  with  a(t)  the  other  can  be  found  from  the  appropriate 
equation  below 


§ 

♦ 


;  2 

-  cos 


2e)1/2 

§2)1/2 


(21) 


(22) 


•  • 

The  appropriate  sign  for  the  derivatives  0(t)  and  <i»(t)  must  be 

prescribed  with  a(t)  since  either  sign  may  produce  feasible  aircraft 

motion.  The  Initial  value  of  the  state  variable  to  be  found  must  be 

defined,  then  numerical  Integration  of  equations  21  or  22  gives  the 

vfi  iue  of  the  state  variable  over  time.  Either  'i'(t)  or  6(t)  can  be 

reiJaced  by  a(t)  in  any  of  the  four  known  sets  of  variables  in  which 

they  appear.  Essentially,  the  new  variable  a(t)  is  added  to  the 

original  11  variables  in  the  problem  and  an  additional  equation, 

either  21  or  22,  is  added  to  the  set  of  equations  to  be  integrated. 

To  solve  equations  14b  and  14c  for  unknown  control  variables  at  any 
•  • 

time,  0(t)  or  <Kt)  are  evaluated  with  the  above  equations. 
Therefore,  the  solution  steps  outlined  to  solve  each  of  the  four 
problems  with  known  <p  and  0  still  apply  when  a(t)  is  p-->8cribed. 
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LIMITATIONS  OF  THE  AIRCRAFT  MODEL 

The  model  assumed  to  govern  aircraft  motion  has  inherent 
limitations  that  must  be  considered  in  its  application.  Specifically, 
values  of  the  prescribed  variables  and  the  resulting  aircraft  motion 
must  be  considered  in  view  of  constraints  on  the  aircraft  motion  not 
implicit  in  the  model  The  moments  an  aircraft  can  physically  achieve 
at  any  flight  condition  are  an  important  set  of  constraints.  Rapid 
changes  in  the  aircraft  attitude  require  large  moments,  so  these  con¬ 
straints  may  be  violated  during  high  angular  rate  aircraft  maneuvers. 
Calculating  aircraft  angular  velocity  during  a  trajectory,  as  given  in 
equation  6,  allows  the  angular  rates  to  be  monitored.  Aircraft  motion 
at  large  angles  of  attack  or  sideslip  generally  cannot  be  predicted 
with  the  trajectory  model  since  problems  in  maintaining  controlled 
flight  can  develop  in  these  flight  regimes.  Additional  constraints, 
due  to  structural  and  engine  operating  limits  can  often  be  represented 
by  a  feasible  flight  envelope  constructed  as  a  function  of  Mach  number 
and  pressure  altitude.  These  are  usually  found  in  individual  aircraft 
operations  manuals.  A  more  detailed  discussion  of  constraints 
frequently  encountered  in  trajectory  analyses  is  presented  in  [2]. 

The  accuracy  of  calculated  fuel  use  will  be  dependent  on  the 
accuracy  of  the  force  and  specific  fuel  consumption  data  as  well  as 
the  accuracy  of  the  approximations  required  for  numerical  solution  of 
the  equations  governing  flight.  Flight  test  results  or  fuel  flow  data 
from  operations  manual  provide  data  to  validate  the  fuel  use 
calculations.  In  certain  applications  of  this  model,  for  example, 


when  comparing  the  fuel  use  of  different  trajectories,  the  relative 
difference  in  fuel  use  values  is  of  primary  importance.  If  absolute 
fuel  flow  values  are  to  be  used,  as  in  mission  planning,  they  should 
be  used  conservatively. 


EXTENSIONS  OF  THE  AIRCRAFT  MODEL 

Equations  14c  through  16  have  been  formulated  to  be  solved  as 
five  different  trajectory  generation  problems.  These  equations 
provide  the  basis  for  solving  several  other  problems  of  interest  in 
aircraft  dynamics.  Flight  path  parameters  or  functions  defined  in 
terms  of  state  and  control  variables  can  be  optimized  using  elements 
of  the  model  described  here.  Typical  problems  include  minimizing  fuel 
flow  with  respect  to  time  or  distance  and  maximizing  climb  rate 
subject  to  a  set  of  constraint  equations.  The  constraint  equations 
may  be  for  example,  the  force  equations  for  symmetric  flight  in  the 
vertical  plane.  A  number  of  optimization  problems  are  formulated  in 
[2].  Example  problems  are  also  formulated  in  [4]  and  numerical 
methods  for  solving  such  problems  sre  presented. 

Problems  in  the  area  of  aircraft  stability  and  control  are 
formulated  in  part  with  the  force  equations  of  motion.  The  moment 
equations  governing  the  angular  motions  of  the  aircraft  are  added  to 
these  equations.  Depending  on  assumptions  about  the  forces  and 
moments,  the  force  and  moment  equations  may  be  coupled  and  require 
simultaneous  solution.  In  the  trajectory  model  presented  here,  the 
equations  are  assumed  independent.  The  moments  acting  on  the  aircraft 


can  be  evaluated  after  a  trajectory  solution  when  validating  the 
feasibility  of  a  high  rate  uaneuver  is  desired*  A  development  of  the 
moment  equations  and  examples  of  their  application  are  presented  in 


IV.  NUMERICAL  METHODS 


The  numerical  algorithms  required  to  solve  equations  lAa  through 
16  have  been  identified  in  the  previous  section.  First,  prescribed 
variables  have  to  be  evaluated  to  solve  the  equations.  In  general, 
these  may  be  arbitrary  functions  constructed  for  a  particular  problem. 
However,  one  method  is  presented  in  this  section  for  simply 
prescribing  variables  as  functions  of  time.  Algorithms  for  solving 
algebraic  equations,  integrating  differential  equations,  and  interpol¬ 
ation  are  also  needed  in  the  solution.  The  algorithms  selected  for 
this  problem  are  described  with  discussion  of  their  application  in  the 
solution  of  equations  14a  through  16.  Finally,  an  expression  for 
approximating  derivatives  of  control  variables  and  the  model  of  the 
atmosphere  to  be  used  are  given. 


PRESCRIBING  KNOWN  VARIABLES 

To  solve  the  equations  of  motion,  four  variables  have  to  be 

defined  as  functions  of  time.  The  variables  that  can  „.e  prescribed 

include  a,  0,  n£,  <j> ,  V,  iji ,  0  and  a.  Not  only  the  values  of  the 

prescribed  variables  are  needed,  but  the  values  •>£  the  time 

derivatives  must  be  calculated  for  all  the  variables  except  a  and 

•  »  • 

nc*  The  derivatives  V,  6,  and  i|>  appear  in  the  acceleration 

•  *  • 

components  of  equations  14a  through  14c  while  a,  0,  and  $  ar.; 
needed  to  calculate  aircraft  angular  velocity  using  equation  6. 


Furthermore,  the  time  derivatives  of  the  prescribed  variables  should 

be  continuous  functions  of  time  to  emure  that  acceleration  and 

angular  velocity  of  the  aircraft  are  continuous.  This  requirement 

arises  from  neglecting  the  aircraft  moment  equations  governing  the 

angular  motions.  Variables  with  linear  continuous  first  derivatives 

can  be  constructed  by  defining  a  sequence  of  constant  second 

derivatives  over  time.  Let  u(t)  denote  a  variable  to  be  prescribed 

•  *• 

with  initial  values  uQ  and  uQ  given  at  time  tQ.  Assume  u^  is  a 

known  sequence  of  time  ordered  second  derivatives  of  u(t)  at 

•  • 

times  tj,  i  -  l,2,...,n.  Rirther  assume  u.  is  constant  on  the 
Interval  t^_,j  <  t  _<  t^ .  Then  u^(t)  and  u^(t)  on  <  t  _<  t^ 

"an  be  found  as 

"  2 
u.At 

u4(t)  -  — ^ — +  ui-lAt  +  ui-l  (23a) 

u±  ( t )  *  ^At  +  u1_1  (23b) 

where  At  ■  t  -  t. .  . 

i.  I 

In  applying  the  above  equations,  u^t)  and  ui(t)  will  represent  the 
value  of  a  control  or  state  variable  and  its  first  time  derivative. 

In  case  of  a(t)  and  nc(t),  the  first  time  derivatives  need  not  be 
continuous,  so  they  can  be  more  simply  prescribed  by  using  equation 
23b  and  equating  u(t)  to  a(t)  or  nc(t).  An  example  of  using  the 
above  method  can  be  illustrated  by  prescribing  a  constant  turn  rate, 


ip.  Suppose  i|»Q  and  are  both  zero  at  tQ.  Then  a  constant  turn 
rate,  ^  could  be  achieved  by  t^,  and  maintained  until  t2  by 
specifying 


Since  u(t)  and  u(t)  on  t^_^  <  t  <  t^  are  evaluated 
Independently  of  t^  in  equations  23a  and  23b,  the  point  in  time 

••  N*  •• 

when  u(t)  switches  from  u^  to  Uj^  does  need  to  be  explicitly 
•  • 

defined  with  u^.  It  will  be  useful  to  allow  t^  to  be  optionally 
defined  as  the  time  when  any  state  variable,  control  variable,  their 
derivatives,  or  acceleration,  say  v(t),  crosses  a  threshold  value, 

n 

c,  during  the  trajectory.  Specifically,  u^  and  u^  are  evaluated 
on  the  interval  <  t  _<_  t ^ ;  t^  -  min(t)  such  that  v(t)  c 

(or  v(t)  c) .  Thus,  instead  of  a  defining  t^  explicitly,  v(t), 
c,  and  the  desired  logical  operator  can  be  defined.  As  an  example  of 

U3ing  this  option,  suppose  an  aircraft,  initially  in  level  flight  at 

velocity  V,  is  to  increase  engine  speed  to  n^,  and  accelerate  to 
V^.  First  initial  conditions  are  defined  for  all  variables  except 
a  and  nc  which  are  solved  at  Lq  by  selecting  the  known  variables 
to  be  u  =•  (V,  <|>,  6,  0),  The  appropriate  initial  conditions  and 

second  derivatives  prescribed  for  steady  flight  maintained  for  one 

. 

second  wili  be  represented  by:  uq  *  (Vq, 0,0,0);  =*  (0,0, 0,0); 

t^  •  1  sec.  The  values  of  the  controls  a  and  nc,  unknown  prior  to 

the  solution,  corresponding  to  the  prescribed  steady  flight  condition 


are  found.  After  switching  to  the  prescribed  variable  set  u  *-  (\Ji , 

0,  n~ ,  2)>  level  acceleration  will  be  accomplished  by  controlling  n 
as  follows: 

u^  =  (0,  0,  n  ,  0);  tj  ■  min(t)  s.t.  n£> 

u3  =  (0,  0,  0,  0);  t^  -  min(t)  s.t.  V(t)  >  V ^ 

••  • 

u4  -  (0,  0,  -nc,  0);  -  min(t)  s.t.  V(t)  <  0 

•• 

u5  =  (0,  0,  0,  0);  t,.  -  5  sec 

•  •• 

Recalling  n£  =*  u(t),  U2  increases  the  engine  speed,  nc,  at  a 

rate  nc-  After  time  t2  the  engine  speed  is  constant  at  a  value 

greater  than  n£.  At  time  t^,  defined  by  V(t)  >_ V,  the  engine 

» 

speed  decreases  at  a  rate  -nc  until  positive  acceleration  ceases, 

• 

V(t)  _<  0.  Then  the  aircraft  maintains  level  flight  for  5  seconds  at  a 
constant  engine  speed.  The  comparison  of  a  variable  with  the 
threshold  value  to  determine  t^  will  occur  at  time  increments  of 
At  equal  to  the  numerical  integration  stepsize.  Therefore,  the  exact 
value  of  a  variable  at  time  t^  cannot  generally  be  predicted  prior 
to  the  trajectory  solution. 

RUNGE-KUTTA  INTEGRATION 

A  fourth  order  Runge-Kutta  algorithm  will,  be  used  to  numerically 
integrate  the  first  order  differential  equations  14a  through  16. 
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Dependent  on  the  problem  formulation,  as  many  as  seven  simultaneous 
equations  must  be  integrated.  Let  y  represent  the  vector  of  state 
variables  to  be  integrated,  x  the  vector  of  prescribed  state 
variables  and  control  variables,  and  f(y,  x)  the  vector  of  ordinary 
differential  equations  in  a  given  problem.  The  control  variables  in 
x  may  include  those  found  as  roots  of  equations  14a,  14b,  and  14c 
so  x  will  be,  in  general,  a  function  of  y  as  w^ll  as  time.  The 
general  integration  problem  can  be  expressed  as 

y  *  f(x,  y),  y(tQ)  -  yQ  (24) 

An  approximation  to  y(t)  at  discrete  points  fci  *  fco  +  iAt» 
i  *  l,2,...,n,  is  desired  where  At  is  a  constant  step  size.  Let 
the  approximation  to  the  solution  y(t^)  be  noted  y^.  The  following 
fourth  order  Runge-Kutta  algorithm  will  be  used  to  calculate  5^+1 • 

yi+i  -  yt  +  +  2k3  +  k4)  (25) 

kj  -  f(yif  x(tj,  yx)) 

*  f(yi  +  j  Atk1,  x(ti+  i  At,  y±  +  y  Ati^)) 

^3  =  ^yi  +TAt^2’  x(ti  +YAt'  yi  +yAt^2^ 

^4  "  ^yi  +  AtS’  x(ti  +  At»  yi  + 


Runge-Kutta  algorithms  are  developed  using  Taylor  series 
expansions  of  tb->.  unknown  solution  y(t).  Neglecting  higher  order 
terras  in  the  expansion  results  in  truncation  error.  An  estimate  of 
the  error,  et,  given  in  [5]  for  fourth  order  integration  is 


e 


t 


16 

15 


<yi+l,2 


yi+l,l> 


(26) 


Here,  ^  is  an  approximation  of  the  solution  y(tj+i)  with 

truncation  error  et  resulting  from  a  step  size  At.  The  term 
y i_l_l  2  is  the  approximation  of  y(ti+i)  based  on  two  integration 
steps  of  size  At/2.  Truncation  errors  can  be  evaluated  at  intervals 
during  integration  and  compared  to  a  threshold  error  value  for  each 
state  variable.  If  the  error  is  exceeded,  the  size  of  At  can  be 
decreased.  The  differential  equations,  f(y,  x),  must  be  evaluated 
four  times  for  integration  across  At.  This  implies  four  evaluations 
of  the  forces  and  specific  fuel  consumption.  Prescribed  variables, 
only  dependent  on  time,  are  to  be  evaluated  at  t^,  t^  +  At/2,  and 
ti'-l*  variables  in  x  which  are  roots  of  equations  must  be  found 

four  times. 


NEWTON-RAP HSGN  ALGORITH-i 

The  Newton-Raphson  algorithm  will  be  used  to  find  the  control 
variables  satisfying  the  set  of  equations  14b  and  14c,  when  tl  j 
velocity  vector  attitude  is  prescribed,  or  14a,  14b,  and  14c  when 
total  velocity  vector  is  prescribed.  In  the  first  case,  the  controls 


(0  ,  a)  or  (<}> ,  a),  must  be  found;  in  the  second  case,  (nc,  0,  a) 
or  (nc>  >  °0  are  the  controls  to  be  found«  Let  f^,  f2  and  £3 
equal  equations  14a,  14b,  and  14c  solved  for  zero  and  x  be  a  vector 
containing  the  unknown  control  variables  required  to  satisfy  f(x)  - 
0.  Given  an  initial  estimate,  Xq,  the  approximation  to  the  solution 
is  iteratively  incremented,  ■  x^  +  Ax^  using 
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f32 

f  33 

f3(X1)_ 

(27) 


where  f -  the  partial  derivative  of  fj  with  respect 

to  control  variable  x^. 


'Jhen  the  velocity  attitude  is  prescribed  using  (\p ,  0,  nc ,  0)  or 
(*K  nc,  4> ) ,  the  problem  reduces  to  solving  two  equations,  f2 

and  f-rj,  for  the  appropriate  controls.  Symmetric  flight  in  the 

\ 

vertical  plane  and  zero  sideslip  flight  also  reduce  the  number  of 
equations  to  be  solved.  When  the  velocity  vector  is  prescribed,  f2 
is  deleted  from  the  problem  for  flight  in  the  vertical  plane  while  for 
zero  sideslip  flight  f2  is  solved  independently  for  $  using 
equation  19.  Similarly,  when  velocity  attitude  is  prescribed,  only 
and  a  appear  in  equation  27  during  symmetric  flight  in  the 
vertical  plane  or  zero  sideslip  flight. 

Multiple  roots  may  exxst  for  any  of  the  sets  of  equations  to  be 
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solved.  Using  different  initial  estimates,  Xq,  when  starting  the 
algorithm  may  provide  the  values  of  the  multiple  roots.  Finding  the 
multiple  roots  may  be  desired  when  starting  the  numerical  solution  at 
steady  flight  conditions,  for  example.  The  Newton-Raphson  algorithm 
requires  that  the  inverse  of  the  partial  derivative  matrix  exist. 
However,  convergence  is  not  guaranteed.  If  problems  arise  other  root 
finding  algorithms  may  be  employed.  Several  alternative  approaches 
are  presented  in  [6].  The  algorithm  stops  when  the  absolute  value  of 
elements  in  the  increment  vector,  Ax^,  are  less  than  the  elements  in 
given  vector  e.  Since  controls  have  to  be  evaluated  four  times 
during  integration  across  At,  the  choice  of  e  will  be  Important  in 
determining  the  time  required  to  integrate  the  remaining  differential 
equations. 

In  the  cases  where  three  equations  are  solved,  nine  of  12 

possible  partial  derivatives  (four  possible  control  variables)  appear 

in  equation  27.  The  expressions  for  the  partial  derivatives  are  not 

presented  here,  but  they  include  the  force  coefficient  functions  and 

partial  derivative  of  these  functions.  B.tth  the  function  values  and 

partial  derivatives  have  to  be  evaluated  by  interpolation  at  each 

iteration.  Specifically,  the  partlals  to  be  evaluated  are  ,  Cq  , 

a  a 

C/i  ,  and  Ki»i  . 
p  “c 

INTERPOLATION  METHOD 

Interpolation  between  discrete  values  of  the  force  coefficient 
and  specific  fuel  consumption  functions  arire  in  the  solution  of 


equations  14a  through  16.  Also,  partial  derivatives  of  these 
functions  with  respect  to  control  variables  will  be  required.  Natural 
cubic  spline  interpolating  functions  will  be  employed  to  meet  these 
requirements.  To  define  these  functions,  assume  n  values  of  a 
function  on  one  dimension,  f(x),  are  given  at  the  base  points  x^, 
i  -  l,2,...,n>  Then  n-1  third  order  cubic  polynomials,  g^(x),  are 
found  by  requiring  continuous  first  and  second  derivatives  on 
(xi,xn).  They  are  uniquely  determined  and  called  natural  cubic 
splines  when  the  second  derivatives  gM(x^)  and  g"(xn)  are  defined 
to  be  zero.  The  remaining  g?  -  g"(x^)  are  found  from  the  solution 
of  n-2  linear  equations  whose  coefficients  are  determined  by  the 
function  values  and  base  points.  The  resulting  coefficient  matrix,  is 
tri-diagonal  and  easily  solved  by  elimination  and  substitution.  The 
problem  formulation  and  solution  methods  are  presented  in  [7]. 

Once  the  values  for  g£  are  found,  then  f(x)  for  x^<x<x^+^, 
is  approximated  as 


.  3  2 

f(x)  "  l(xi+rx)  ’Axi(xi+rx)] 


«"<W 


6Ax,  (<x-xi>3  -  AxJ(x-Xl)] 


f(*i)  f(x,  .) 

+  it-  <xi+rx)  +  tt-  (x-v  • 


(26) 


where  -  xi+l_xi* 


The  derivative,  f’(x),  follows  from  this  equation. 
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To  apply  the  cubic  spline  interpolation  to  functions  defined  on 
two  dimensions,  cubic  splines  will  be  calculated  along  both 
coordinates  for  every  base  point.  A  function  f(x^,x2)  requires  two 
second  partial  derivatives,  3  8ij/3x£»  k  “  1*2,  at  each  base  point 
(xi  X2  j).  These  need  only  be  calculated  once  and  stored  with  the 
original  function  values  for  each  bast  cvln  Equation  28  can  be  used 
to  approx  mate  an  arbitrary  point  f(xj_,  where  (xj_,  x2)  is  In 

the  rectangle  £.  xi  <_  j+- ,  i  —  *‘i  l*.  x2  j+l»  a  cubic 

spline  parallel  to  one  coordinate,  say  x2,  is  defined.  First  the 
two  points  f(x^,  x2>j)  and  X2,J+1)  are  found  with  equation 

28  applied  twice  along  edges  of  the  rectangle  parallel  to  the  x^ 
coordinate  using  the  stored  function  values  and  partial  derivatives 
with  respect  to  x^  Next,  two  second  partial  derivatives  in  the 
x2  direction  3  g(xlt  x?  j)/3x2  and  BgCx^  x2,j+i)/3x2  are  needed* 
These  will  b'i  approximated  by  linear  interpolation  between  partial 
derivatives  in  the  x2  direction  known  at  the  corner  points  of  the 
rectangles . 
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(X1,1+1  "  Xl,i) 


k  -  j,  j+1 


Once  these  second  derivatives  are  calculated,  then  the  two  function 
values  and  the  two  derivatives  required  to  evaluate  f(x^,  x2)  using 


pressure  and  density  are  found  as  functions  of  altitude  by  solving  the 
hydrostatic  equation  and  the  ideal  gas  law  simultaneously.  The  speed 
of  sound,  a,  can  be  accurately  modeled,  see  [8],  as  a  function  of 
absolute  temperature,  r,  using  reference  values  Sq  and  xQ  by 

a  -  a0(t/T0)1/2  (31) 

A  standard  atmosphere  is  frequently  used  in  modeling  aircraft 
flight.  The  NACA  standard  atmosphere  as  presented  in  [8]  will  be 
approximated  by  assuming  a  sea  level  temperature  of  59°F  decreasing  at 
a  rate  of  3.56  x  10"3oF/ft  to  -67.6°F  at  about  35.3  kft,  where  the 
temperature  remains  constant  to  well  above  conventional  jet  aircraft 
ceilings.  At  59°F,  the  sea  level  sound  velocity  is  1,117  ft/sec,  and 
the  assumed  pressure  is  2116.2  lb/ft  *  With  this  information,  the 
properties  of  air  can  be  found  and  tabulated,  then  linear 
interpolation  will  be  used  for  quickly  evaluating  the  air  properties 


V.  SUMMARY  AND  CONCLUSIONS 


To  assess  aircraft  tactics  using  quantitative  models  of  weapons 
likely  to  engage  an  aircraft,  a  quantitative  description  of  an 
aircraft's  trajectory  is  required.  In  this  report,  a  model  conven¬ 
tionally  used  in  aircraft  trajectory  analyses  is  formulated  as  a 
system  of  seven  equations  and  eleven  variables.  The  variables 
governed  by  the  equations  define  the  aircraft's  position,  velocity, 
fuel  use,  and  attitude.  A  particular  aircraft's  capabilities  are 
represented  by  the  force  and  specific  fuel  consumption  functions  that 
appear  in  the  equations.  It  is  assumed  these  functions  can  be  modeled 
by  dimensionless  coefficients  defined  on  two  dimensions. 

To  determine  a  unique  trajectory  governed  by  the  equations,  a  set 
of  variables  must  be  prescribed  as  functions  of  time.  Different  sets 
of  variables  can  be  selected  to  conveniently  prescribe  different 
segments  of  a  trajectory  depending  upon  the  aircraft  flight  condition 
or  maneuver  desired.  Either  the  control  variables,  the  attitude  of 
the  velocity  vector,  or  the  attitude  and  magnitude  of  the  velocity 
vector  can  be  prescribed  with  the  latt-L  two  sets  also  including 
selected  controls.  The  derivatives  of  the  velocity  attitude  angles 
determine  the  flight  path  norr,ial  acceleration,  a  quantity  useful  for 
prescribing  maneuvers-  One  of  the  velocity  attitude  angles  can  be 
replaced  by  normal  acceleration  in  the  set  of  prescribed  variables  if 
an  additional  equation  governing  the  replaced  variable  is  added  to  the 
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model.  To  ensure  aircraft  acceleration  and  angular  velocity  are 
continuous,  prescribed  state  and  control  variables  should  possess  at 
least  continuous  first  derivatives.  Evaluating  prescribed  variables 
using  a  sequence  of  constant  second  time  derivatives  is  one  possible 
procedure  that  meets  the  continuity  requirements. 

Numerical  solution  of  the  equations  involves  application  of 
algorithms  for  integrating  first  order  differential  equations,  finding 
roots  of  nonlinear  equations,  and  interpolation  to  approximate 
function  values.  During  integration  across  a  time  interval,  both  the 
root-finding  and  interpolation  algorithms  must  be  applied  at  each 
derivative  evaluation.  Integration  error  can  be  estimated  and  the 
time  interval  reduced  if  desired.  Each  of  the  algorithms  will  be  best 
implemented  independently  of  the  equations  to  be  solved  and  ir. 
separate  subroutines.  The  logic  required  to  solve  the  system  of 
equations  will  be  reflected  in  the  sequence  of  subroutine  calls  and 
the  parameters  that  are  passed  at  each  call.  This  approach  will 
reduce  the  effort  necessary  to  solve  the  equations  with  different  sets 
of  aeroo/namic  force  and  engine  data.  Adding  arbitrary  functions  for 
prescribing  variables  cr  extending  the  methods  to  solve  other  dynamics 
problems  is  simplified  if  the  algorithms  and  equations  of  the 
trajectory  generation  model  are  well  defined  in  the  software 
Implementation. 

Aircraft  trajectories  found  using  the  model  in  this  report  must 
be  evaluated  with  respect  to  external  constraints  on  the  feasibility 
of  the  aircraft  motion.  These  include  limits  on  the  aircraft 
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structure,  engine  operation,  and  aircraft  control.  The  operations 
manuals  for  a  particular  aircraft  provide  information  on  the  flight 
limitations  that  must  be  observed.  If  maintaining  aircraft  control 
during  a  flight  path  is  questioned,  the  implied  moments  can  be. 
calculated  and  compared  to  the  maximum  control  moments  available.  The 
aircraft  moment  coefficient  data  and  moment  of  inertia  properties  /, re 
needed  for  these  calculations. 
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